function [R_tr,R_re,R]=corre_matrix(Tr,Rr,cor)
[R_tr,y]=createR(cor(1),Tr);
[R_re,y]=createR(cor(2),Rr);
R=kron(R_re,R_tr);
R=sqrtm(R);
